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ABSTRACT 

We study the effect of extragalactic point sources on satellite observations of the cosmic mi- 
crowave background (CMB). In order to separate the contributions due to different foreground 
components, a maximum-entropy method is applied to simulated observations by the Planck 
Surveyor satellite. In addition to point sources, the simulations include emission from the 
CMB and the kinetic and thermal Sunyaev-Zel'dovich (SZ) effects from galaxy clusters, as 
well as Galactic dust, free-free and synchrotron emission. We find that the main input com- 
ponents are faithfully recovered and, in particular, that the quality of the CMB reconstruction 
is only slightly reduced by the presence of point sources. In addition, we find that it is pos- 
sible to recover accurate point source catalogues at each of the Planck Surveyor observing 
frequencies. 

Key words: methods: data analysis - techniques: image processing - cosmic microwave 
background. 



1 INTRODUCTION 

A new generation of cosmic microwave background (CMB) satel- 
lite missions are currently in the final stages of design. The NASA 
MAP satellite is expected to be launched by 2000, followed by the 
ESA Planck Surveyor in 2007 (Bersanelli et al. 1996). Both mis- 
sions will provide detailed all-sky maps of the CMB anisotropies, 
leading to definitive measurements of the CMB power spectrum. 
This should allow tight constraints to be placed on fundamental 
cosmological parameters and distinguish between competing the- 
ories of structure formation in the early Universe such as inflation 
and topological defects. 

The maps produced by these satellites will, however, contain 
contributions from various foreground components, most notably 
Galactic dust, free-free and synchrotron emission as well as the 
kinetic and thermal SZ effects from galaxy clusters. In addition, 
significant contamination from extragalactic point sources is also 
expected. It is therefore clear that in order to obtain maps of the 
CMB anisotropies alone, it is necessary to separate the emission 
due to these various components. 

In a previous paper, Hobson et al. (1998) (hereafter Paper I) 
use a non-linear maximum-entropy method (MEM) to separate the 
emission due to the different foreground components from simu- 
lated Planck Surveyor observations of a 10 x 10 deg 2 field (see 



Bouchet et al. 1997; Gispert & Bouchet 1997). It was found that 
faithful reconstructions may be produced not only of the CMB 
anisotropies but also of the Galactic components and the kinetic and 
thermal SZ effects from massive clusters. The analysis presented in 
Paper I does not, however, include the effect of extragalactic point 
sources, since no reliable simulations were available at that time. 

In this paper, we apply a generalised version of the MEM sep- 
aration algorithm to simulated Planck Surveyor observations that 
also include a population of point sources. This allows us to as- 
sess the effect of point sources on the accuracy with which the 
CMB anisotropies and foreground components can be recovered. 
The simulated maps of the point sources were created for each 
Planck Surveyor observing frequency using the updated number 
counts in each frequency band predicted by Toffolatti et al. (1998). 
These predictions are based on the Danese et al. (1987) model for 
the evolution of radio selected sources (adopting an average spec- 
tral index a = for compact flat-spectrum sources up to ~200 
GHz and a break to a = 0.7 at higher frequencies) and on the 
model C of Franceschini et al. (1994) for the evolution of far-IR 
selected sources, updated by Burigana et al. (1997) to account for 
the isotropic sub-mm component estimated by Puget et al. (1996). 
The point source map at each frequency was produced by Poisson 
distributing all the sources in the 10~ 5 Jy < 5(v) < 10 ly flux range. 
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Unless sources are identified and subtracted down to faint flux lim- 
its, the contribution of clustered sources to the CMB anisotropics is 
generally found to be small in comparison with the Poisson term. 
As discussed by Toffolatti et al. (1998), far-IR selected sources 
should dominate the number counts of bright objects at frequen- 
cies v > 300 GHz, whereas radio selected sources should dominate 
at lower frequencies. The two underlying parent populations are 
largely different and uncorrelated, consisting of normal spirals and 
starburst galaxies in the far-IR and mostly compact AGN, blazars 
and quasars in the radio. 



2 THE SEPARATION ALGORITHM 

Our aim is to reconstruct the CMB anisotropics and the foreground 
components in the presence of both instrumental noise and con- 
tamination due to extragalactic point sources. A straightforward 
initial step might be to subtract the point-source contribution from 
those pixels that are obviously contaminated. The procedure can, 
however, lead to difficulties, since a particular contaminated pixel 
can usually only be identified in a few of the frequency channels. 
Therefore, since the spectral behaviour of point sources is generally 
rather complicated, any interpolation between frequencies must be 
performed with caution. An alternative Fourier-filter approach to 
removing point sources from CMB maps has been proposed by 
Tegmark & de Oliveira-Costa (1998). 

In our approach, we do not remove any contaminated pixels 
from the simulated data. Instead, we consider the point sources as 
an additional contribution to the 'noise' and allow the generalised 
version of the MEM algorithm to reconstruct the CMB and fore- 
ground emission in the presence of the point sources. The MEM 
technique is discussed in detail in Paper I and so here we will out- 
line the basic points and instead concentrate on how to address the 
effect of point sources within the existing formalism. 

At any given frequency v, the total rms temperature fluctua- 
tion on the sky in a direction x is given by the superposition of the 
physical components. In Fig. [l] we show the 10 x 10 deg 2 simu- 
lated input maps for all the different physical components (except 
the point sources) at a reference frequency Vo = 300 GHz; each 
map consists of 400 x 400 pixels and has been convolved with a 
Gaussian beam with a FWHM of 4.5 arcmin, which is the highest 
angular resolution of the Planck Surveyor satellite. The cell size for 
each map is 1.5 arcmin. 

If we observe the microwave sky in a given direction x at iif 
different frequencies, we obtain an n y-component data vector that 
contains the observed temperature fluctuations in this direction at 
each observing frequency plus instrumental noise. The vth compo- 
nent of the data vector in the direction x may be written as (see 
Paper I) 

d v{x) = Y,Pv{\x-xj\) £ F vp s p (xj) +r\ v (x) +£ v (x), (1) 
j=\ P =i 

where N p denotes the number of pixels in each of the simulated 
input maps shown in Fig. In this expression we have chosen to 
separate the contributions from the six physical components shown 
in Fig. |j| and the point sources. The former are collected together 
in a signal vector with n c = 6 components, such that s p (x) is the 
signal from the pth physical component at the reference frequency 
Vo = 300 GHz. The corresponding total emission at the observing 
frequency v is then obtained by multiplying the signal vector by the 
>ifXri c frequency response matrix F vp that includes the spectral be- 
haviour of the considered components as well as the transmission 



of the vth frequency channel. This contribution is then convolved 
with the beam profile P v (x) of the relevant channel. Since the spec- 
tral dependence of the point sources is very complicated, we cannot 
factorize their contribution in this way and so they are added into 
the formalism an extra 'noise' term. Thus r) v is the contribution 
from point sources as observed by the Planck Surveyor satellite at 
the frequency v. Finally, e v is the expected level of instrumental 
noise in the vth frequency channel and is assumed to be Gaussian 
and isotropic. Figs g and ^ show the simulated maps at each ob- 
serving frequency for the LFI and HFI respectively. The simulated 
data are identical to those presented in Paper I, but also include the 
point source simulations of Toffolatti et al. (1998). 

The assumption of a spatially-invariant beam profile in (|l]) al- 
lows us to perform the reconstruction more effectively by work- 
ing in Fourier space, since we may consider each fc-mode indepen- 
dently (see Paper I). Thus, in matrix notation, at each mode we 
have 



d = Rs + J7 + e = Rs + C, 



(2) 



where d, rj and e are column vectors each containing n y complex 
components and s is a column vector containing n c complex com- 
ponents. In the second equality we have combined the instrumental 
noise vector e and the point-source contribution rj into a single 
vector 'noise' vector £. The response matrix R and has dimensions 
iif xn c and its elements are given by R vp (k) = P v (k)F vp . 

As discussed in Paper I, the MEM formalism should not itself 
induce correlations between elements of the reconstructed signal 
vector. However, the elements of the signal vector s at each Fourier 
mode may well be correlated, this correlation being described by 
the n c x n c signal covariance matrix C defined by 



C(k) = (s(k)s'(k)), 



(3) 



where the dagger denotes the Hermitian conjugate. At each Fourier 
mode, the pth diagonal element contains the value of the ensemble- 
averaged power spectrum of the pth physical component at the ref- 
erence frequency Vq. Moreover, if prior information is available 
concerning these correlations, we would wish to include it in our 
analysis. We therefore introduce the vector of 'hidden' variables h, 
related to the signal vector by 



s = Lh. 



(4) 



where the n c x n c lower triangular matrix L is obtained by per- 
forming a Cholesky decomposition of the signal covariance matrix 
C = LL T . The reconstruction is then performed entirely in terms of 
h and the corresponding reconstructed signal vector is subsequently 
found using (Q). 

Using Bayes' theorem, we choose as our estimator h of the 
hidden vector to be that which maximises the posterior probability 
given by 



Pr(h|d)cxPr(d|h)Pr(h) 



(5) 



where Pr(d|h) is the likelihood of obtaining the data given a par- 
ticular hidden vector and Pr(h) is the prior probability that codifies 
our expectations about the hidden vector before acquiring any data. 

We assume an entropic prior probability for the hidden vector 
h of the form 



Pr(h) oc exp[<xS(h,m) 



(6) 



where S(h, m) is the cross entropy of the complex vectors h and m, 
where m is a model vector to which h defaults in absence of data. 
The form of the cross entropy for complex images and the Bayesian 
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Figure 1. The 10 x 10 deg 2 realizations of the six input components used to make simulated Planck Surveyor observations: (a) primary CMB fluctuations; (b) 
kinetic SZ effect; (c) thermal SZ effect; (d) Galactic dust; (e) Galactic free-free; (f) Galactic synchrotron emission. Each component is plotted at 300 GHz and 
has been convolved with a Gaussian beam of FWHM equal to 4.5 arcmin, the maximum angular resolution proposed for the Planck Surveyor. The map units 
are equivalent thermodynamic temperature in /jK. 
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Figure 2. The 10 x 10 deg 2 maps observed at each of the four observing frequencies of the Planck Surveyor Low Frequency Instrument (LFI): (a) 30 GHz, (b) 
44 GHz, (c) 70 GHz, (d) 100 GHz. At each frequency we assume a Gaussian beam with the appropriate FWHM and a sampling rate of FWHM/2.4. Isotropic 
noise with the relevant rms has been added to each map. The map units are equivalent thermodynamic temperature in /jK. 



method for fixing the regularising parameter a are discussed in Pa- 
per I. 

Let us now turn our attention to the form of the likelihood 
function Pr(d|h) in Bayes' theorem In the presence of instru- 
mental noise alone, it is reasonable to assume that the noise con- 
tribution is well described by a Gaussian distribution, but in our 
case the noise term also contains a contribution from the point 
sources, and so we might expect that a Poisson distribution would 
be more appropriate for this component. However, since we are 
performing the reconstruction in the Fourier domain, the noise on 
each fc-mode will contain contributions from a wide range of scales. 
Therefore, provided the point sources are distributed on the sky in a 
statistically-homogeneous manner, we would expect from the cen- 
tral limit theorem that the noise in the Fourier domain is described, 
at least approximately, by a Gaussian distribution. Therefore, we 
assume that the likelihood function is given by 

Pr(d|h) oc exp^-C + N _I c) 

oc exp[-(d-RLh) + N _I (d-RLh)] (7) 

where in the last line we have used (Q). The noise covariance matrix 
N has dimensions rif x n/ and at any given A>mode is given by 

N(k) = (C(k)C(k)). (8) 

Therefore, at a given Fourier mode, the vth diagonal element of 
N contains the ensembled-averaged power spectrum at that mode 
of the instrumental noise plus the point source contribution to 
the vth frequency channel. The off-diagonal terms give the cross- 
correlations between different channels; if the noise is uncorre- 
lated between channels, only the point sources contribute to the 
off-diagonal elements. 

The argument of the exponential in the likelihood function (^) 
may be identified as (minus) the standard % 2 misfit statistic, so we 
may write Pr(d|h) oc exp[— % 2 (h)]. Substituting this expression, to- 



gether with that for the prior probability given in (|6|), into Bayes' 
theorem, we find that maximising the posterior probability Pr(h|d) 
with respect to h is equivalent to minimising the function 

3>(h)=x 2 (h)-aS(h,m). 

This minimisation can be performed using a variable metric min- 
imiser (Press et al. 1994) and requires only a few minutes of CPU 
time on a Sparc Ultra workstation. 



3 APPLICATION TO SIMULATED OBSERVATIONS 

We now apply the MEM analysis outlined above to the simulated 
Planck Surveyor data shown in Figs ^ and |3| Clearly, the quality 
of the reconstructions will depend on our prior knowledge of the 
input components. Following Paper I, we assume that the frequency 
spectrum behaviour of the components is accurately known, but we 
consider two extreme cases relating to our knowledge of the power 
spectra of the input components and the confusion noise due to 
point sources. 

For the six input components shown in Fig. Ill any prior 
power spectrum information may be incorporated into the algo- 
rithm through the signal covariance matrix C given in (^). Strictly 
speaking, since we are in fact reconstructing the vector of hidden 
variables h, rather than the signal vector s, this information actually 
resides in the matrix L given in (Q). Similarly, any knowledge of 
the correlation structure of the instrumental noise and point source 
contamination can be included via the noise covariance matrix N. 

3.1 Reconstructions with full power spectrum information 

As our first case, we assume knowledge of the azimuthally- 
averaged power spectra of all six input components in Fig. to- 
gether with the azimuthally-averaged cross power spectra between 
components; these contain cross-correlation information in Fourier 
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Table 1. The rms of the residuals per 4.5 arcmin FWHM Gaussian beam 
(in fiK) of the MEM reconstructions in the presence and absence of point 
sources, assuming full power spectrum information. 



Component 


e rms no sources 


firms with sources 


CMB 


5.90 


7.50 


Kinetic SZ 


0.85 


0.86 


Thermal SZ 


3.90 


4.40 


Dust 


1.60 


3.80 


Free-Free 


0.30 


0.38 


Synchrotron 


0.05 


0.06 



space, so that the matrix L is fully specified. We also assume that 
the noise covariance matrix N is fully specified as follows. In addi- 
tion to the (diagonal) part of the noise matrix due to the instrumen- 
tal noise, we calculate an empirical estimate of the contribution (to 
the diagonal and off-diagonal elements) of the point sources. This 
estimate is calculated directly from 10 independent realizations of 
the point source maps at each of the 10 Planck Surveyor frequency 
channels. 

The resulting MEM reconstructions of the physical compo- 
nents are shown in Fig. ^ and are plotted on the same greyscale 
as in Fig. [l] to allow a straightforward comparison with the true 
input maps. We see that the main input components are faithfully 
reconstructed, except for a few contaminated pixels containing the 
brightest point sources. Perhaps most notable is the fact that the 
majority of contaminating point sources present in the data maps, 
shown in Fig. ^ and are not present in the reconstructions. 

In particular, the CMB has been reproduced extremely accu- 
rately and is virtually indistinguishable from the true input map, 
apart from a spurious bright point source in the lower left-hand 
corner of the reconstruction. As we might expect the dust emission 
is also accurately recovered, and contains none of the numerous 
point sources present in the high-frequency data maps shown in 
Fig. ^. The free-free emission, which is highly correlated with the 
dust, has also been reconstructed reasonably accurately and dis- 
plays no evidence of point source contamination. The recovery of 
the synchrotron emission is, however, somewhat poorer. Although 
the main features of the true input map are recovered, the recon- 
struction contains obvious residual point sources. Perhaps most im- 
pressive is the reconstruction of the thermal SZ effect, since the 
MEM algorithm has been successful in reproducing the effect in 
most of the bright clusters and has misidentified only a few point 
sources as clusters. At the reference frequency of 300 GHz, these 
misidentified point sources appear mostly as negative features. The 
kinetic SZ effect has also been reconstructed in a number of clusters 
having bright thermal SZ effects. In fact, the quality of the kinetic 
SZ reconstruction is similar to that obtained in Paper I. 

In order to quantify the effect of point sources on the accuracy 
of the reconstructions, in Table |l| we compare the rms residuals, 
fiims, per 4.5 arcmin FWHM Gaussian beam, of the reconstructions 
in Fig. Q with those obtained in Paper I in the absence of point 
sources. As we might expect, for each component the value of e rms 
is slightly larger when the effects of point sources are included, but 
the errors are still reasonable small. In particular, we note that the 
rms error on the CMB reconstruction only increases from 5.9 to 7.5 

As an intermediate case, reconstructions were also performed 
assuming that only the diagonal elements of the noise covariance 
matrix were non-zero. This therefore includes information about 



the power spectra of the instrumental noise and point source maps 
at each observing frequency, but ignores any cross-correlation of 
the point source between frequencies. In this case, the quality of the 
reconstructed maps is only slightly worse and, in particular, the rms 
of the residuals in the CMB reconstruction increases marginally to 
8.3 ^K. 

In Fig. ^| we plot the power spectra of the reconstructed maps 
(faint lines) and compare them with the power spectra of the true 
input maps (bold lines) as shown in Fig. [jj. We see from Fig. || that, 
in spite of the contamination due to point sources, the power spec- 
trum of the MEM reconstruction of the CMB component (panel a) 
is accurate up to £ ~ 2000, at which point the reconstruction be- 
gins to underestimate the true power spectrum. For the kinetic SZ 
effect (panel b) the power spectrum of the reconstruction underes- 
timates the true spectrum at all multipoles, but for the thermal SZ 
effect (panel c) the reconstructed spectrum follows the true one up 
to I ~1000. The most accurately reconstructed power spectrum is 
that of the Galactic dust emission (panel d), which closely follows 
the true spectrum up to I ~ 3000. Finally, the power spectra of the 
free-free and synchrotron reconstructions (panels e and f) slightly 
underestimate the true spectra across all multipoles. It was further 
found that assuming a diagonal noise covariance matrix did not ap- 
preciably alter the results. As discussed in Paper I, the errors on 
the reconstructed power spectra may be easily calculated, and we 
find that in each case the 68 percent confidence limits on the recon- 
structed power spectra encompass the true power spectrum at all 
multipoles. 

3.2 Reconstructions with no power spectrum information 

For our second case, we take the opposite view to that adopted 
above and assume that almost no power spectrum information is 
available. This corresponds to assuming a flat (white-noise) power 
spectrum for each component out to the highest measured Fourier 
mode. The levels of the flat power spectra are chosen so that the to- 
tal power in each component is approximately that observed in the 
input maps in Fig. [lj Furthermore, the noise covariance matrix con- 
tains only the diagonal contribution form the instrumental noise, 
and no attempt is made to include the effects due to point source 
contamination. As discussed in Paper I, when no power spectrum 
information is assumed, the MEM algorithm is iterated by using 
the current reconstructions to update the ICF matrix L before em- 
barking on the next iteration. The solution was found to converge 
after 8 such iterations. 

The resulting reconstructions are shown in Fig. ^| In this case, 
it is most encouraging that the CMB reconstruction is still very 
accurate, even for this very pessimistic scenario. Once again the re- 
construction is, at least by eye, almost indistinguishable from the 
true map shown in Fig. [I], except for the presence of a very bright 
point source in the lower left-hand comer of the map. We also see 
that the Galactic dust component is accurately recovered and con- 
tains very few of the point sources present in the high-frequency 
channel maps shown in Fig. ^. The thermal SZ effect has also 
been recovered, but only in the brightest clusters, and the recon- 
strucion also contains numerous misidentified point sources that 
appear as sharp negative features. The reconstructions of the re- 
maining three components are very poor. In particular, the kinetic 
SZ reconstruction has simply defaulted to zero in the absence of 
any useful data, whereas the reconstructions of the Galactic free- 
free and synchrotron emission are heavily contaminated with point 
sources and no obvious features in the true input maps have been 
recovered. 
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Figure 4. MEM reconstruction of the 10 x 10 deg 2 maps of the six input components shown in Fig. [l|, using full power spectrum information (see text). The 
components are: (a) primary CMB fluctuations; (b) kinetic SZ effect; (c) thermal SZ effect; (d) Galactic dust; (e) Galactic free-free; (f) Galactic synchrotron 
emission. Each component is plotted at 300 GHz and has been convolved with a Gaussian beam of FWHM equal to 4.5 arcmin. The map units are equivalent 
thermodynamic temperature in /jK. 
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Table 2. The rms of the residuals per 4.5 arcmin FWHM Gaussian beam 
(in /jK) of the MEM reconstructions in the presence and absence of point 
sources, assuming nopower spectrum information. 



Component 


£rms no sources 


e lms with sources 


CMB 


6.10 


8.70 


Kinetic SZ 


2.68 


2.70 


Thermal SZ 


4.35 


5.80 


Dust 


1.90 


4.80 


Free-Free 


0.44 


1.60 


Synchrotron 


0.07 


0.23 



In Table g we compare the rms residuals, erm S , per 4.5 arcmin 
FWHM Gaussian beam, of the reconstructions shown in Fig. ^ with 
the corresponding results obtained in Paper I in the absence of point 
sources. We see that the value of the e lms has increased for all the 
components with respect to the case when point sources are not 
present in the simulations. Nevertheless, the relative errors for the 
CMB, thermal SZ and Galactic dust components are still reason- 
ably small. In particular, we note that even in this pessimistic sce- 
nario, the rms error on the CMB reconstruction is still only 8.7 /liK. 

Finally, in Fig. [7] we plot the power spectra of the reconstructed 
maps and compare them with the input power spectra. For the CMB 
(panel a), we see that the reconstructed power spectrum again fol- 
lows the true spectrum very closely up to t ~ 1800, at which point 
it drops directly to zero. The dust power spectrum (panel d) is also 
very accurately recovered up to I ~ 2500. For the thermal SZ ef- 



fect (panel c) follows the true spectrum up to £ « 200, but then 
begins gradually to overestimate the true spectrum before drop- 
ping to zero at I ~ 1200. The additional power in the reconstructed 
power spectrum is a result of misidentifying point sources as ther- 
mal SZ effects from clusters. As expected from the reconstructed 
maps shown in Fig. ^| the power spectra of the true kinetic SZ, 
free-free and synchrotron emission are not recovered and the recon- 
structed free-free and synchrotron power spectra oscillate widely 
about the true spectrum before dropping to zero at I ~ 300. We 
again find that the true power spectrum in each case lies within the 
estimated one-sigma errors on the reconstructed power spectrum at 
all multipoles. 

We note that for the CMB and dust components the power 
spectra are apparently recovered over a larger range of multipoles 
than when point sources were not included in the analysis (see Pa- 
per I). The reason for this becomes clear if we include full power 
spectrum information for the six input components shown in Fig. [lj, 
but do not include the point source contribution in the noise co- 
variance matrix. In this case, the reconstructed power spectra of 
the CMB and dust are overestimated and possess a distict bump 
in the multipole range t 1500-2500 and £ « 2500-3500 respec- 
tively. This feature is due to the contribution to high multipoles of 
the point sources. We must be careful when interpreting the power 
spectra of the recovered maps at high multipoles, since these an- 
gular scales may be contaminated by point source emission despite 
the fact that only a few point sources are obvious in the recon- 
structed maps. 
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Figure 6. MEM reconstruction of the 10 x 10 deg 2 maps of the input components shown in Fig. |lj using no power spectrum information (see text). The 
components are: (a) primary CMB fluctuations; (b) kinetic SZ effect (c) thermal SZ effect; (d) Galactic dust; (e) Galactic free-free; (f) Galactic synchrotron 
emission. Each component is plotted at 300 GHz and has been convolved with a Gaussian beam of FWHM equal to 4.5 arcmin. The map units are equivalent 
thermodynamic temperature in ,uK. 
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100 1000 « 100 1000 - 1 100 1000 
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Figure 7. The power spectra of the input maps (bold line) compared to to the power spectra of the maps reconstructed using MEM with no power spectrum 
information (thin line). In panel (b) only the input power spectrum is shown for the kinetic SZ effect, since the reconstructed map in this case is everywhere 
zero. 



4 RECOVERY OF POINT SOURCE CATALOGUES 

We have so far concerned ourselves primarily with the reconstruc- 
tion of the CMB and other foreground emission in the presence of 
point sources. Nevertheless, an additional aim of the Planck Sur- 
veyor mission is to compile point source catalogues over a wide 
range of frequencies. In this section we therefore investigate how 
accurately we can recover the input point source population used in 
the simulations. 

It is clear that the generalised MEM approach discussed in 
Section ^ has in no way been optimised to achieve the goal of re- 
covering point sources. Indeed, the mechanism for preventing the 
inclusion of sources into the reconstructed CMB and diffuse fore- 
ground maps was to regard the sources as an additional noise com- 
ponent. Nevertheless, given the absence of point sources in the re- 
constructions shown in Figs ^ it is clear that some information 
concerning these sources can be recovered by comparing the recon- 
structions with the input data. 

We use the following straightforward technique to obtain in- 
formation about the point source distribution at each Planck Sur- 
veyor observing frequency. We first make simulated data at each 
frequency using the MEM reconstructions as inputs (as opposed 
to true maps shown in Fig. [I]). For each observing frequency the 
MEM reconstruction of each physical components is projected in 
frequency using its known spectral behaviour. The resulting maps 
are then added together to produce the MEM reconstruction of the 
sky at this frequency in the absence of point sources. This map 
is then convolved with the appropriate Planck Surveyor observing 
beam to produce mock data at this frequency. These mock data are 
then subtracted from the true input data (which contain the point 



source contribution). Since the MEM reconstructions contain very 
few point sources, we thus hope to recover an approximate point 
source map at each Planck Surveyor observing frequency. 

Using the MEM reconstructions shown in Fig ^, which assume 
full power spectrum information, we find that the recovered point 
source maps are quite accurate. For the sake of brevity we display 
only the reconstructions for the 44 and 353 GHz channels, which 
are shown in Figs ^ & |^ respectively. These channels are chosen 
simply because they lie near the centres of the frequency ranges 
of the LFI and HFI respectively. We emphasise that these maps 
represent the point sources as observed by the Planck Surveyor, i.e. 
convolved with the appropriate observing beam. 

In Fig. [| we see that at 44 GHz the reconstructed point source 
map in panel (b) is very similar to the input map in shown in panel 
(a). Indeed all of the prominent point sources have been recovered 
at approximately the correct flux. There is, however, some evidence 
of residual diffuse emission in the top left-hand corner of the recov- 
ered map. We quantify the point source recovery by identifying and 
extracting point sources from the input and recovered maps using 
the SExtractor algorithm (Bertin & Arnouts 1996). This algorithm 
first fits an unresolved background to an image and then identifies 
point sources superposed on this background. Using the standard 
default settings, SExtractor fitted an unresolved background with 
mean level 0.05 Jy per beam area to both the true and recovered 
point source maps shown in Fig.|^(a) and (b). After subtracting this 
background, the resulting rms fluctuation was found to be 0.01 Jy 
for each map and SExtractor then identified sources down to this 
level. Thus the detection flux limit in each case was 0.06 Jy. For the 
input map, 93 sources were found as compared to 85 for the recov- 
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Figure 8. (a) The input point source map at 44 GHz convolved with the ap- 
propriate Planck Surveyor observing beam, (b) The recovered point source 
map at 44 GHz. (c) A histogram of point source fluxes for the input map 
(heavy line) and recovered map (faint line) obtained using SExtractor; see 
text for details. The greyscales in (a) and (b) are flux density in Jy 



ered map. The resulting histograms of source counts for the input 
map (heavy line) and the recovered map (faint line) are shown in 
Fig. pic) and are clearly very similar. 

In Fig. ^| we show the results for the 353 GHz. Since the 
FWHM of the observing beam at this frequency is only 4.5 arcmin, 
the number of distinguishable point sources is much higher than at 
44 GHz. We note once more that the input and reconstructed point 
source maps are very similar. Indeed, even very faint sources are 
accurately recovered. SExtractor again fitted an unresolved back- 
ground component with mean value 0.05 Jy per beam area to both 
maps, resulting in a residual rms fluctuation of 0.01 Jy in each case. 
Thus, as for the 44 GHz channel, the detection flux limit was 0.06 
Jy. The algorithm identified 1869 sources in the true map and 1676 
sources in the recovered map. The histograms of the point source 
counts are shown in Fig. ^(c) and are seen to agree reasonably well. 
There is, however, some evidence that the recovered map has too 
few sources in the range 0.10-0.12 Jy and too many in the range 
0.06-0.08 Jy. This is most likely caused by some of the flux as- 
sociated with the point sources being incorrectly assigned to very 
low-level diffuse foreground emission by the MEM separation al- 
gorithm. 

The results at other observing frequencies are of similar qual- 
ity to those shown in Figs |s| & []. We thus find that the input point 
source catalogues at each Planck Surveyor frequency can be re- 
covered to reasonable accuracy using the above method. We re- 



Figure 9. As for Fig. |8| but for the 353 GHz Planck Surveyor frequency 
channel. 

emphasise that this has been achieved with no direct subtraction of 
point sources from the data, but instead by the application of the 
MEM separation algorithm in which the point sources were con- 
sidered as an extra 'noise' contribution. 

The recovered point source maps discussed above were cre- 
ated using the MEM reconstructions obtained by assuming full 
power spectrum information. Nevertheless, we find that the results 
are not appreciably affected by using the MEM reconstructions ob- 
tained assuming no power spectrum information shown in Fig. ^. 
At the lower LFI frequencies the flux of the brightest sources is un- 
derestimated by about a factor of two since some of flux has been 
misidentified by MEM as synchrotron or free-free emission. In the 
HFI frequency range, however, the dominant diffuse emission is 
due to the CMB and Galactic dust, which are both accurately re- 
constructed, and the resulting point source recovery is of a similar 
quality to that shown in Fig. EH 



5 DISCUSSION AND CONCLUSIONS 

In this paper we study the effect of point sources on the reconstruc- 
tion of CMB anisotropies and foreground emission from simulated 
satellite observations. In particular, we apply a generalised form of 
the maximum-entropy method (MEM) developed in Hobson et al. 
(1998) to analyse simulated observations by the Planck Surveyor 
satellite of a 10 x 10 deg 2 patch of sky. 

We do not attempt to remove contaminated pixel from the 
data, but instead introduce information about the power spectrum 
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of the point sources in each frequency channel and the correlations 
between frequencies. The resulting reconstructions show that, even 
in the presence of contaminating point sources, the CMB emission 
can still be recovered with an accuracy of about 8 fjK, even for 
the pessimistic scenario in which we assume no prior information 
about the covariance structure of the sky emission. It is possible that 
this accuracy could be improved by a suitable scheme for removing 
point sources from the data before the MEM component separation 
algorithm is applied. Moreover, given some prior knowledge of the 
power spectra of the CMB and foregrounds, not only can the accu- 
racy of the CMB reconstruction be improved, but it is possible to 
obtain accurate reconstructions of the thermal SZ effect from clus- 
ters and the Galactic dust, free-free and synchrotron emission. The 
kinetic SZ effect may also be recovered in clusters that have a large 
thermal effect. 

As a final point, we find that by using the MEM reconstruc- 
tions of the CMB and foreground components it is possible to sub- 
tract simulated mock data from the true data in order to recover the 
point source contribution at each of the Planck Surveyor observing 
frequencies. The recovered point source maps can then be analysed 
using, for example, the SExtractor algorithm in order to compile 
point source catalogues. 
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